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o 

' We calculate the exact degeneracy of states corresponding to the area operator in 

(N 

, the framework of semiclassical loop quantum gravity, using techniques of combina- 

D ■ 

pLn I torial theory. The degeneracy counting is used to find entropy of apparent horizons 



(N 



derived from generalised coherent states which include a sum over graphs. The 



^ . correction to the entropy is determined as exponentially decreasing in area. 

o 
o 

o 

^ : I. INTRODUCTION 

o 

O^' One of the most enigmatic and interesting physical quantities is the entropy of a horizon. 

As evident now, the horizons may or may not enclose singularities. From the discovery of 



the irreversibility associated with horizon area in any classical process, and laws of black 
hole mechanics [1-3], entropy of a horizon has become synonymous with thermodynamic 
entropy. This identification requires the search for a theory of quantum gravity and a 
corresponding 'microscopic' counting of entropy. In the absence of a complete theory of 
quantum gravity, the horizon area has been arbitrarily quantised in units of Planck length 
squared and a degeneracy of 2 assigned to every bit of area [2]. Later, reduced phase 
space quantisations actually achieved an 'area quantisation' [4], and the area spectrum 
was identified as equispaced numbers, in terms of the fundamental Planck length squared. 
However, these quantisations were not complete 'quantum theories' and inclusion of the 
full degrees of freedom of gravity might have completely different conclusions from reduced 
phase space quantisations. Thus, the search for an explanation for the microscopic origin 
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of what is now known as Bekenstein-Hawking entropy, which precisely equals A///4/p 
continues. This search for the origin of entropy is quite obviously linked to the search for 
a complete theory of quantum gravity. Within the realities of the four-dimensional world, 
the most promising theories are the standard techniques of path-integral quantisation and 
canonical quantisation of gravity. String theory, the more radical approach, remains in the 
realm of the unreal, far from a quantum description of a 4-dimensional non-supersymmetric, 
non-charged black hole [5]. In the path-integral approach, due to difficulties in actual 
evaluation of the functional integral, one has to discretise the space-time. This discretisation 
of space-time for a Schwarzschild black hole in the dynamical triangulation framework 
appears in a recent work [6]. However, it is only in the framework of canonical quantisation 
of gravity that a derivation of entropy, including the resolution of singularity at the center 
of black hole, has been obtained, mainly prompted by the discovery of a quantisation 
of area. For reduced phase space quantisations using geometrodynamics, there has been 
recent progress in [7] . For reduced phase space resolution of black hole singularities in loop 
quantum gravity, see [8]. However, as shown in [9], inclusion of full degrees of freedom of 
quantum gravity might change the conclusions completely. Thus, quantisation of horizons 
is still a subject of research. Since we are working in the loop quantum gravity framework, 
much of the formalism is motivated by a previous derivation of entropy by counting the 
states of a Chern Simons theory induced on the horizon boundary [10, 11]. However, the 
entire programme as described in [11] of entropy is characterised by the following 

a) The microstates for entropy are associated with only the horizon, and the rest of the 
space-time does not manifest itself in any way. 

b) In actual black hole processes the horizon is not a boundary of space-time, and informa- 
tion is known to flow across the horizon in a quantum mechanical description. 

c) The area eigenstates at the horizon are taken as exact eigenstates, which is a good 
approximation in the large area, limit, but should not be the correct physics, if the horizon 
exists only classicaly. 

d) The actual counting after the imposition of the horizon boundary conditions on the state 
were not exact, and the Immirzi parameter remains a mystery [12, 13]. 

However, the calculation was remarkable in the sense that the quanta of area of the 
horizon was given by 4,ti^ j{j + 1)/^ with the j corresponding to a spin degree of freedom. 
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For j = 1/2 this would precisely correspond to the 2 degrees of freedom associated with unit 
area of 1^. This would make the previous derivations by physicists a reality. The calculation 
also showed the way for later calculations of black hole entropy in the loop quantum gravity 
framework. 

One of the recent developments which has a very correct perspective is to represent an 
entire black hole space-time by a semi-classical or coherent state [14, 15]. Since black holes 
are essentially classical solutions of Einstein's equations, and the light cone is a classical 
concept, horizon entropy should be associated with semi-classical or coherent states. The 
coherent state is defined as a function of complexified phase space variables, and the expec- 
tation value of operators are closest to their classical value in this state [16, 17]. The horizon 
is realised as a solution to a difference equation, which is exact only in the classical limit 
[18]. The entire wavefunction for the coherent state is defined for a spatial slice of the black 
hole which includes the horizon as well as the singularity. The singularity is resolved in the 
coherent state, and is recovered only in the classical limit. Most importantly, the entropy can 
be calculated from first principles as Sbh = — Trplnp, where Sbh is the entropy associated 
with the horizon, and p is the density matrix. This density matrix was derived by tracing 
over the coherent state wave-function inside the horizon. Thus the origin of entropy is sim- 
ilar to an entanglement entropy calculation, first observed in [19]. Thus it is not surprising 
that the entropy is proportional to the area of the surface which bounds the set of states 
which are integrated out. What is completely new here is that this entanglement entropy 
is completely 'gravitational' in origin and hence finite. Thus previously infinite calculations 
of entanglement entropy for scalar fields in black hole space-times [20] are now completely 
finite and correct. 

However, the coherent state discussed in [16] calculated the entropy for one particular 
graph. Here, we discuss the situation where, within certain symmetry restrictions, there is 
a sum over inequivalent graphs. The kinematic Hilbert space is a tensor sum of individual 
Hilbert spaces [21]. Each sub-Hilbert space has a 'coherent state' corresponding to that 
graph. This interpretation is necessary, as within the formalism of [15], the coherent state 
is defined over graph dependent variables. The 'generalised' coherent state is defined in the 
tensor sum of Hilbert spaces, and the operators are also tensor sum of operators for each 
graph. For the semi-classical limit, one has to do a suitable averaging process. 

To obtain entropy, the density matrix is obtained for each graph coherent state separately. 
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The density matrix is a tensor sum over density matrices for each graph. In the final 
expression, the entropy is the logarithm of the degeneracy counting associated with the 
degeneracy of the area operator at the horizon. This counting is done using generating 
function techniques of combinatorics. The final answers are remarkably simple, and exact. 
The main observations of our paper are 

i) Entropy is proportional to of the horizon. 

ii) The final value of the Immirzi parameter is given as 8 ln(3/2 + 75/2) = 8 ln(2.618). This 
value is the same value as was obtained in [22], where the equispaced area spectrum used 
here was used to find the entropy. The ensemble of spins was taken to be made up of 
completely distinguishable elements. 

iii) The corrections to horizon entropy are actually exponentially decreasing with area, in 
contrast to previously observed logarithm of area corrections [23]. The log area correction 
appears only if the cardinal number of edges inducing the horizon taken to be fixed. 
However ab-initio, there is no reason to restrict the number of spins to a particular value. 

Having said all these, one must emphasise that the entropy derivation is done using co- 
herent states defined on the kinematical Hilbert space, as opposed to the physical Hilbert 
space where the diffeomorphism and the Hamiltonian constraints are satisfied [24, 25]. This 
would be a very interesting formalism to develop. What we hope is that the generating 
function techniques of combinatorics used here to calculate the degeneracy of states and the 
'semi-classical' realisation of entropy using coherent states will be useful for future calcula- 
tions of entropy. The derivation of the exponentially decreasing correction terms to entropy 
is a very new result, and it will be very interesting to determine if this is the theoretical 
constraint which will get the correct theory of quantum gravity. 

In the next section we give the formalism for the derivation of generalised coherent states, 
and the sum over Hilbert spaces interpretation of the classical limit. The third section 
contains the derivation of entropy from a exact combinatorial counting, and finally the 
fourth section ends with conclusions. 

II. GENERALISED COHERENT STATES 

In the new- variables or Sen-Ashtekar-Barbero-Immirzi variables, the gravitational canon- 
ical variables of the spatial three metric and the corresponding momentum (function of ex- 
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trinsic curvature) are re- written in terms of densitised triads and a SU(2) gauge connection. 
The defining equation for these variables are 

(3E^E'j6'' = detq q'^' = - (3K', (1) 

where Ej are the densitised triads, qab the three metric, the SU(2) connection. is 
the spin connection corresponding to the triads, and = —7= — e'^^Kah has the information 
about the extrinsic curvature Kab associated to the spatial slicing. (3 is the Immirzi pa- 
rameter, and represents a one parameter ambiguity in the definition of the variables. This 
ambiguity creates inequivalent quantisations, and hence the value determined with the black 
hole entropy is in some sense a experimental result for the theory. The actual quantisation is 
done in the variables of holonomy and momentum defined on graphs comprised of piecewise 
analytic edges, and their corresponding dual graph, comprised of 2-dimensional surfaces. 
The basic phase space variables are now: 

K{A) = Pexp( J A.dx) Pi = ^Tr[T^/ie ^ hpEh'^h^^] (2) 

where /le(^) denotes the holonomy along one edge e and is a path ordered exponential of 
the gauge connection along the edge. P/ is the corresponding momentum evaluated on a 
dual graph comprised of 2-surfaces and is a function of the densitised triads as well as the 
holonomy of the gauge connection along the edge e which cuts the dual surface, and the 
hp denote the holonomy along the edges defined on the dual 2-surface S. a is a dimensional 
constant, which is determined from the length scales of the theory. 

Thus in the quantisation of the above variables, 'graphs' have a very important role. 
In generic quantisation, the Hilbert space is taken as defined over the configuration space, 
which is the projective limit of graphs. This is possible as the graphs form a directed set 
7' > 7 where 7' can be obtained by sub-dividing the edges of 7. The projection operators 
Py^ then map the graph configuration spaces. The kinematic Hilbert space is defined over 
the space of the connections modulo the gauge transformations. The states in this space 
are not diffeomorphism invariant, and neither do they satisfy the Hamiltonian constraint. 
However, they are relevant, as they have a correct semi-classical limit, and the expectation 
value of operators are closest to their classical values, evaluated on a classical metric. 

Now, these variables have a Poisson algebra which can be lifted to commutators as a 
quantisation. The quantised Hilbert space is given by orthonormal basis states Ijmn > 
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for each edge. A complexification of the phase space takes the variables to SL(2,C) valued 
elements 

= exp-^^'^'/' K (3) 

(T^ are the generators of SU(2) and are ^ where a are Pauli matrices) and the Hilbert 
space wavefunction can be defined in terms of these in the Segal-Bergmann representation. 
The kernel of the transformation, which is also called a coherent state transform, is the 
coherent state itself, as originally defined by Hall. This of course is a function of both 
heiQei where now ge is interpreted as a classical label for the phase space point at which the 
wavefunction is peaked. Now, given a graph F comprised of edges, the complete coherent 
state for the entire graph is given by a tensor product over coherent states for each edge. 
Thus, 

^*(^?e,M =n^*(^?e,M (4) 

e 

{t = Ip/a is the semiclassical parameter.) Imposing the trivial representation of intertwiners 
at the vertices would give a gauge invariant state, and the expectation values of the oper- 
ators can be determined using the same techniques as above, and they are also closest to 
their classical values [26]. The coherent states written down in [16-18], were defined for a 
Schwarzschild black hole, also in a particular choice of coordinates, and one particular choice 
of graphs. The entropy was also defined for one particular graph. 

However, there is no reason to choose one particular graph: the graphs can be varied. 
Two graphs, which belong to a directed set, can be said to be equivalent, as the expectation 
values of the holonomy and the momenta can be obtained by taking the product of the 
holonomy and the momenta of the other graph. In a directed set 7' > 7, the edges of 7' 
are obtained by subdividing the edges of 7. Since, for a edge e, which is broken into two 
edges 61,62, the holonomy would be = he^he^, one can always obtain operators of 7, if one 
knows the operators of 7'. However, the reverse is not true, and hence, given a projective 
family one can always identify a 'minimal graph'. The kinematic Hilbert space is then a 
tensor sum of Hilbert spaces for these minimal graphs. 

Each such Hilbert space can be characterised by a coherent state, and a corresponding 
classical limit obtained for the operators as defined in terms of the expectation values in 
those states. The generalised coherent state would be defined on the tensor sum of Hilbert 
spaces. However, there is a subtlety about the graphs and the classical limit of operators 
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defined on them, which first has to be completely understood. The graphs in consideration 
for the classical limit have to be regarded in terms of two separate questions. 
l)The actual choice of a 'minimal' arbitrary graph might be classified as follows. 

1) The number of vertices. 

ii) The number of edges linked at a particular vertex, or the valency of the graph. 

iii) The orientation of the edges. 

iv) The total number of edges, which can in principle be countably infinite. 

v) Any symmetries of the graph. An example of this would be a cubic graph, which is 
six valent, and can be very useful in determining 3-dimensional semi-classical limits, as 
considered in [27]. 

2) The actual embedding of the graph in a classical metric. These would be decided by the 
following considerations: 

i) The classical metric. 

ii) The edge lengths of the graph as measured in the classical metric. 

iii) Symmetries of the classical metric. 

iv) The holonomy and the momentum, which will characterise the physical phase space. 

For example, a cubic graph can have a spherical embedding, as well as a non-spherical 
embedding. Thus, there is no unique way of embedding a given graph in a classical space 
time, and hence the same 'quantum graph' can correspond to different 'classical graphs'. 
In the coherent states for gravity, the actual classical graph can be thought as labelled by 
the Qe (3). Thus, the semi-classical limit corresponding to a given Hilbert space will be 
characterised by the graph, as well as the set of labels denoted as {ge}- 

The graph as discussed in [16] for spherically symmetric classical space-times was actually 
fixed by the symmetries of the classical metric. The graph was chosen so that the nth vertex 
coincided with the 0th vertex. Further, each vertex was 6-valent. While taking the 'classical 
limit' of this graph, what remained fixed was the cardinal number of edges, and vertices, 
and the above stated 'spherical symmetry'. Thus, each coherent state, corresponding to a 
given graph, will now be characterised by 

^r,{,.}' (5) 

where F denotes a minimal graph. 
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In case of the spherically symmetric embedding, the choice of edges along the coordinate 
lines of r, 6*, (p of the sphere ensures that the graph samples maximal points near the origin. 
However, here, the number of vertices is fixed per volume, and hence the density decreases 
as one approaches the asymptotics. To sample more points in the exterior region, new 
non-radial edges have to be introduced. The questions of singularity resolution and entropy 
calculation are however quite relevant as given in [18]. Also, a very similar choice of graph 
seems to be sufficient to describe the interior of a reduced Schwarzschild space-time in [8]. 
Thus, even with a rather careful choice of graph, there remains the ambiguity of the classical 
graph. 

The question is, given a graph F, how does one evaluate the sum over classical labels 
(yfg? For coherent states, there is a very natural measure for the classical labels, given by the 
Liouville measure, and this was discussed in [18]. The answer obtained there is that the area 
induced by the edges given by ^fWP^ = Pe has to be very small. This can be interpreted 
as an indication of a natural way to ensure microscopic sampling. However, in a generic 
situation, a graph can induce all possible values of Qe- In the very particular restrictive case 
we consider here, we give each graph a unique classical set of labels {Qs}- If we take the 
Liouville measure seriously, then, as considered in [18], only those graphs which induce the 
minimum area of the 2-surfaces would be preferred. In the coherent states area spectrum, 
this would correspond to j = or semiclassical area t. This would not yield any entropy 
associated with the horizon. However, given any graph, we choose the classical graph to 
sample the minimal areas induced. 

Thus a generalised 'coherent state' will be: 

^ = ®r,9.^r,c;e (6) 

The symbol ® denotes cartesian product, and the operators will be tensor sums. The 
eigenvalues will then be just sums of eigenvalues for each space. Now, we concentrate on the 
'annihilation operators' for such Hilbert spaces, and determine their action on the coherent 
states. This annihilation operator is defined as 

= e''/'e-'^^^^/% (7) 

In the tensor sum Hilbert space, this would correspond to having only the block belonging 
to this particular graph edge as operator valued in the tensor summed operator. Hence, this 
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operator of course is specific to a particular edge, and hence would give a non-zero value, 
iff a particular graph containes the edge. The local properties would be very specific to the 
particular graph, and hence the edges. 

^e|^>=^e|^> (8) 

Thus, only the particular graph edge corresponding to that particular operator will give a 
non-zero eigenvalue. The variable can then be used to estimate the values of the holonomy 
and the momenta using the following: 

e-'"^'"' = gegle-''/' (9) 

and 

= e-^'l^e^P'^'l-'g, (10) 

Now in the generalised coherent state, one has the opportunity of improving the estimate 
obtained from a single graph, which might not sample all possible points of the classical met- 
ric. If one takes the annihilation operator which is the tensor sum of annihilation operators 
of all the sub-Hilbert spaces, then in a given region i?, at least one of the edges will sample 
the points. Thus for the 'global' annihilation operator, which will be a tensor product of 
annihilation operators for edges of each individual graph, there will be a maximal sampling 
of the classical points in the tensor summed Hilbert space. 

There will, however, be a overcounting for 'macroscopic' operators, whenever the same 
edge samples the same classical points in more than one graph. 

This will be evident, in the 'global operator' or the tensor product of the operators for 
each complete graph. 

Let us explain what we mean by a 'global operator'. The annihilation operator corre- 
sponding to a entire graph would be 

= ®ge- (11) 

This operator will have eigenvalues which are 2x2 matrices for each edge. Thus if the 
graph has n edges, the tensor product will be a 2r2 x In matrix. 

If there corresponds another minimal graph, which has m edges, the corresponding 'global 
operator' would be a 2m x 2m matrix. 
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The tensor sum of the operators would be to simply add the 2n x 2n and the 2m x 2m 
matrices to obtain: 

gr,2n® gr,2m- (12) 

Now, what is mainly interesting is the situation where the tensor sum of operators is 
being used to measure classical quantities. If the two graphs have very similar embeddings 
in the classical metric, then since many of the edges sample common classical points, the 
eigenvalues will be repeated. For measurement of macroscopic quantities, e.g. area of a 
2-sphere of radius R, both the graphs would give the same result to first order in the semi- 
classical parameter. 

In the situation considered in [18], this can be the area of the horizon. Each of the graphs 
will give the 

J2e Pe = ^5 and they should all be equal. Thus, though the number of 
edges will differ from graph to graph, due to the macroscopic area restriction, the total area 
has to be same. Thus for the total area operator defined for the 2 graphs above, in the limit 
t 0, 

<-^\A®A\%l)>=2A (13) 

Thus, the operator has to be divided by the number of graphs in the tensor sum to obtain the 
correct value. The point we are trying to make is that the local properties can be recovered 
from any set of edges, and their expectation value on a coherent state adapted to that 
graph. However for operators, charting macroscopic or global regions of the manifold, the 
tensor sum of operators will sample the same classical region more than once, and hence one 
has to find a suitable averaging process. Now, in principle, there can be countably infinite 
graphs, and hence the averaging procedure could be a 'infinite renormalisation' process, 
which physicists are quite familiar with. However, as we discuss, the symmetry restrictions 
of the classical metric might restrict the number of graphs, and thus there will be a method 
to count the graphs. 

Now we specialise to the case of the spherically symmetric apparent horizon considered 
in [18], and observe what a generalised coherent state entropy would correspond to. 

The horizon area is induced by radial edges crossing the horizon, as was motivated from 
the apparent horizon equation in [18]. This is really a very restrictive example of inducing 
area of one particular surface by a set of edges. The number of degrees of freedom which 
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can vary is really the number of radial edges, and their spacing on the surface. 

Let us represent that by a planar circle, and see how the graphs can vary. There are 
actually different ways to change the graph. These can be of the following form 
(i) The cardinal number of edges changes. Here, as explained in the following diagram, 
the introduction of two new radial edges (in bold and in color) crossing the horizon, is 
actually changing the previous minimal graph to another one, as clearly the new graph is 
not obtained by subdividing previous edges. Hence the two Hilbert spaces corresponding to 
the graphs would have to be summed over. 




Horizon 



(ii)The graphs have the same number of edges, however these are not symmetrically 
embedded. Hence the spin assignments jd will differ. In other words the classical labels 
Qe would be different. This would however still be a 'different' minimal graph, which we 
explain subsequently. 




Horizon 
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The reason we are concentrating on the graph at the horizon is because entropy is only 
a function of the edges at the horizon, and these are a countable set. Hence all graphs 
which have the same number of edges with the same set of classical labels Qe at the horizon, 
can be grouped together as entropically 'equivalent class' of graphs and these will be quite 
definitely countable. Hence for the purposes of the calculation of black hole entropy, we will 
sum over only those graphs which differ at the horizon. 

To understand the tensor product structure using a scheme given in [28], we find that 
the Hilbert space can be written as a tensor sum of Hilbert spaces for each graph a, with 
a restriction that H'^ consists of cylindrical functions, which are orthogonal to every other 
graph strictly contained in a. Thus H'^ is a tensor sum over states with non-zero j for each 
edge, and non-zero representation /' at each vertex. Now, let us look at the following graphs 
at the horizon, and identify the 'minimal graphs'. Let us start with a graph with n radial 
edges crossing the horizon. If one tries to reduce the number of edges crossing the horizon, 
since the embedded graph is six-valent, one has to remove an entire set of vertices along 
one radial coordinate. If the graph has different valences away from the horizon, removal 
of vertices immediately outside the horizon, and the vertices within would have the same 
effect. Now, let us start from a graph with n edges, and with a set of removal or addition 
of vertices, see how different graphs can be obtained from the same. For example, removal 
of the ith. edge, i — 1th edge, i — 2th edge gives a graph of n — 3 edges, however, this graph 
would be inequivalent to the graph obtained by removing the ith edge, i — 4th edge and 
i — 8th edge from the same graph. Graph I would differ from Graph II in the distribution 
of the areas induced by the edges at the horizon 2-surface in the classical limit. One thing 
is very clear, since the edges are embedded in the classical metric, the areas induced by the 
edges would correspond to 'distinguishable' spins. The Hilbert space of graph I would have 
the same number of edges as graph II crossing the horizon, however, they are inequivalent 
graphs and their semi-classical limits are different. 
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i-4 




i-8 V 
GRAPH II 



Now, one would think that the density matrix will be a tensor sum over the density 
matrices for each graph. The tracing process of edges 'inside' the horizon is graph specific. 
Hence, each Hilbert space will be associated with its own density matrix, and the combined 
density matrix will be a tensor sum of the density matrices for each graph 

P = ®Pr.,u.y (14) 

Now, thinking of this as derived from a superposition of coherent states for different graphs, 
this is equivalent to ignoring interference effects between two different graph Hilbert spaces. 
We calculate the entropy corresponding to the above, and find the leading term to be 
Bekenstein-Hawking. The corrections are exponentially decreasing in area. However, if one 
includes those, one gets a larger entropy. The leading term is again Bekenstein-Hawking, 
and the corrections are exponentially decreasing in area. We now give a much more detailed 
description of the derivation of entropy. 
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III. ENTROPY OF HORIZONS 

The entropy associated with spherically symmetric apparent horizons has been derived 
in [18]. The apparent horizon equation shows that the 'classical correlations' across the 
horizon have information only about the total areas induced at the horizon by the edges. 
The individual components remain undetermined. This is the origin of degeneracy in 
the density matrix. Each individual Hilbert space has its own density matrix, however, the 
total entropy should include a sum over all possible graphs. As obtained in [18], the density 
matrix for each graph is of the form: 

Pr,{34 = H Pe (15) 

i=l..n 

where the edges are precisely the edges crossing the horizon and inducing the same with 
tiny bits of area, given in terms of the spin as Ph = (jd + where t = l^/a. Each 
pe is diagonal, with the entries being l/(2jci + 1)- We give the derivation here, for the 
density matrix for one particular graph, which has say n radial edges crossing the horizon. 
Ab initio, the tensor product of the coherent states peaked at the individual edges does 
not have information about the continuity of the classical metric across the vertices. This 
has to be included in the definition of the semi-classical limit of the quantum theory, as 
the classical metric is derived from a local differential equation. Hence, the metric of two 
neighbouring points x'^ and + e are obtained by solving the equation of motion, or the 
differential equation corresponding to the Einstein Lagrangian. Thus the holonomy of edge 
ei is related to the holonomy of edge 62, if they are linked at the same vertex. However, as 
in [18], we are interested in determining the correlations in the wavefunction which exist due 
to the existence of an apparent horizon equation. The classical metric beyond the horizon, 
and before the horizon, do not contribute in any way to the correlations across the horizon. 
Since the trace is taken over the wavefunction inside the horizon, the correlations which 
matter are the ones which link the wavefunction of the edges outside the horizon with those 
immediately inside the horizon. This is simply obtained by identifying the apparent horizon 
in the spherically symmetric coordinates, and then identifying the correlations by lifting the 
equation to an operator equation in he, P/, which is also a discretised version of the same. 
Thus the classical metric of space-time is chosen in a particular set of coordinates where the 
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spatial or constant time slices are spherically symmetric and have intrinsic curvature: 



c?s3 = dr^ + r\de^ + sin 



(16) 



The entire curvature of the space-time metric is contained in the extrinsic curvature or Kab 
tensor. Now, if there exists an apparent horizon somewhere in the above spatial slice, then 
that is located clS db solution to the equation 



VaS" + KabS^S^ -K = 



(17) 



where S"' , {{a,b = 1,2,3) denote the spatial indices) is the normal to the horizon. If the 
horizon is chosen to be the 2-sphere, then S"' = (1, 0, 0), and the apparent horizon equation 
function of the metric reduces to: 



Krril 



q' 



-K. 



Keeq' 



- n 

i 0^ — u 



(18) 



The first term of the equation disappears trivially as 1 = q^'^ for any point in the spatial 
slice. Hence the extrinsic curvature along the radial direction plays absolutely no role in the 
apparent horizon equation. In the discrete version of the above equation, the information 
about the extrinsic curvature is contained in the holonomy along the radial edge. The h^,^ 
play absolutely no role in the correlations induced in the coherent state wavefunction due 
to the apparent horizon equation. The radial wavefunction is thus independent and is not 
restricted by the apparent horizon equation. The rest of the equation is converted to a 
difference equation in the variables h^g and /ig^, Peg and Pe^ by a suitable regularisation 
along edges. This equation, as described in detail in [18], is given by 



PL P. 



I r 



V 



Tr 
1 d 



T'h-'[K„v] 



VP op 



Tr[T' ^he] ^PL 



Tr 



Tr 
1 d 



T'h-}[K,,V'I'] 
Tr[T^ ^Kg] PP^^ 



Tr 



T'h:l[he,V'/'] Tr[T^K,[K„V'/'] 

J f 2 J 

0) = (19) 



where quantities evaluated at a vertex Vi outside the horizon get related to quantities inside 
the horizon at a vertex V2. The appropriate regularisations which have been used here are 
given by 

(20) 



e{4>)m) 



ee(0)^e(</,)> 



p^ 
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ei^,)=TrrKl,{he,,,,,V]] (21) 

and 

Kk^) = P^Tr[T^Ko] (22) 

It is quite difficult to explicitly solve (19), but one can gather some information from them, 
in the following way. Symbolically (19) can be written as 

A[B\vi) - B\v2)]C\vi) -D = (23) 
B\v^)C'{v^) - ^ = B\v,)C'{v,) (24) 

Now, though the equation is completely independent of the direction / in the internal gauge 
space, the individual components of B\v2) are dependent on the components of B\vi) as 
projected along the vector . Thus, the apparent horizon equation solves even the internal 
gauge components of Qe or the group valued classical label of the edges ending at vertex V2 
in terms of the classical labels at vertex Vi. 

Now, the apparent horizon equation (19) is thus a difference equation relating the canon- 
ical variables of holonomy and momenta at vertices Vi with vertices f 2 which are angular 
in nature. Inspired by this, and as sketched before, the graph at the horizon is particularly 
chosen to be comprised of only radial edges which cross the horizon, and link vertices of the 
type vi with vertices of the type V2- When the above is lifted to an operator equation, and 
the expectation value of the operator equation evaluated in the coherent states, one ends up 
correlating the wavefunction outside the horizon with those inside. In the t — > limit, this 
is manifested as a relation of the classical labels outside with those inside. 

{9evJ (25) 

This correlation is introduced in the coherent state wavefunction by a conditional proba- 
bility function /((yfg viiOe v2\Ph), which takes care of the classical limit only, and is the same 
as a number in the t limit. This is done to retain the form of the coherent state 
wavefunction as defined in [24], and still have information about the correlation. For the 
calculation of the density matrix, the coherent state can be written as the tensor product 
of Hilbert spaces associated with graph edges which end at a vertex vi outside the horizon. 
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the edges which start at a vertex inside the horizon V2, and the radial edges en which hnk 
the vertices vi and V2- The wavefunction, in the momentum representation, can be written 
in the area eigenbasis which are \jmn > where j is the SU(2) Casimir quantum number and 
m,n = —j... + j. 

I*>= E i^{jo}JH,{ji}\{jo} > \{Jh} > > (26) 
{jo},iH,{ii} 

where the {jo}, {ji} denote all the combined Casimirs of edges outside and edges inside, 
and jn is the horizon Casimir eigenvalue. Now the density matrix formed from them can 
be written as 

p=|$><^| (27) 

If the coefficients 4'{jo},jH,{ji} factorise into products of coefficients, then any reduced density 
matrix obtained from them by tracing over one of the Hilbert spaces will be again a pure 
density matrix. However, if it does not factorise, then the resultant density matrix is a 
mixed density matrix. In the formulation of the coherent state in [24], these would have 
factorised. However, now, due to the apparent horizon equation, these do not factorise, and 
the reduced density matrix is a mixed density matrix. Now, in the t — > limit, the density 
matrix simplifies to just a diagonal one, with only non-zero elements in a (2jc/ + l) x [2jci + l) 
block, where jd is the 'Casimir' corresponding to the classical value of the area induced by 
the radial edge on the horizon. In the t ^ limit, the wavefunction can be written in the 
form 

I* >= E i'{jo}i9{o})i^{jH}i9H)i'{ji}{gi)f{jo},{jH},{ji}(9o,PH,gi)\{jo} > \{jh} > \{ji} > 

{jo},{jH},{jl} 

(28) 

Where the wavefunction is in a semifactorised form, with the information about the corre- 
lation only in the term /, which prevents a complete factorisation of the coefficients. If the 
Qi, go is any different from that determined by the apparent horizon equation of (19), then 
the function is vanishing, otherwise it is non-vanishing. Now, the main point to note here 
is that since the apparent horizon equation is independent of the holonomy of the edges 
at the horizon, and depends only on the eigenvalue Ph, the / function does not have any 
information of the eigenvalues rriHjnH associated with the horizon edges. However, this is 
clearly not true for internal quantum numbers m^o},^{o}y^{i}y^{i} as these determine the 
internal component along the gauge directions, and hence by (24) have to be determined. 
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Thus the degeneracy of the density matrix in the t — > hmit arises only from the horizon 
edge quantum numbers. 

Now, using the exphcit form of the coefficients ^jmn = G~^''^''~^^^^'^T^j{g)mn, and the fact 
that these are non-zero in the classical limit only for certain values of {jo}, {Jh}, {ji}, the 
density matrix has diagonal terms only for 

Here the mH,nH = —jfj-.-jfj. Thus, to obtain Trp = 1, |/p = 2p+T- This completes the 
simple derivation. Thus, in the end, the entropy derived from the above is the number of 
ways to induce the horizon area, in the semi-classical limit. One has just given a framework 
for obtaining the same, using a coherent state wavefunction, and a method of tracing the 
wavefunction within the horizon. The points to note here are, of course, the coherent states 
are still kinematical states, and the Hamiltonian constraints have not been imposed. This 
work still requires much more effort, and work is in progress. The main observation here is 
that entanglement entropy can be completely gravitational in origin, and when so, can be 
completely finite. 

Now the above is the density matrix obtained locally, or due to the tracing over the 
wavefunction for one particular radial edge. The complete density matrix would be given 
by (15). 

The entropy function for the density matrix would be 

Sbh = -Tr(plnp) (30) 

Now, if the number of spins were i = l..n, then the density matrices would be n(2ji + 
1) X n(2ji + 1) (dropping the classical labels) matrices with diagonal entries l/n(2ji + 1), 
and hence the entropy would be of the form: 

Sbh = EM2j^ + l) 

i 

= ln(n(2j. + l)). (31) 

i 

In addition there will be the constraint that the spins sum up to the total area of the horizon 
+ 1/2)^ = ^n/^p- Similarly when the 'generalised' density matrix is considered, the 
density matrix is not obtained as a tensor sum of the above density matrices for each 
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individual graph. That would be a mixed sum, as the interference terms of the different 
Hilbert spaces would not be there. However, we just consider the same here, 

p = Opr (32) 

The combined density matrix has block diagonal entries each of dimension Y\i(2ji + 1) x 
YliC^ji + 1) dimensions with diagonal elements of l/Hii'^ji + !)• The trace of the density 
matrix has to be normalised to be 1. This implies the diagonal elements still have to be 
divided by 1/d, where d is the total number of configurations. The entropy is then 

S = -Tr(p In p) = ^ Eli^(n(2ji + 1)) + Ind (33) 

Ji ^ 

The leading contribution to entropy appears from the Ind term, where d is the total number 
of configurations. This counting can be exactly done and equals 2^"^, where = 2/1/^//^. 
The other term can be shown to be bounded by a constant independant of Ah- This gives 
the entropy to be linear in area. Both of these calculations are carried out in (III A) below. 

To see what a more 'democratic' density matrix might look like, we now derive the density 
matrix from an ab-initio 'generalised' coherent state, which is a state in the tensor sum of 
Hilbert spaces. 

The coherent state for a generic graph is a tensor product of coherent states for each edge 
comprising the graph. The 'generalised coherent state' is then a tensor sum of the above 
coherent states. 

I^'gen >= ©rl^r > (34) 
The density matrix written for the above would be as given below: 

Pgen l^gen ^gen|- (35) 

The above does not break up into a tensor sum of the following 

Pgcn = ©rl^r >< ^rl (36) 

but there remain cross terms in the metric, which can give rise to transition amplitudes 
from one graph to another. This is difficult to physically justify at the semi-classical level, 
however, at any finite order in t where there will be quantum corrections, such transitions 
are completely possible. A given graph can have an edge added to it due to an operator 
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acting on it; for example, the Hamiltonian operator adds edges to a given graph. Thus 
to zeroth order in t, the density matrix should be completely diagonal and as t becomes 
finite, it will have off-diagonal terms also appearing. Now, to calculate the terms, the basic 
procedure would be the same as outlined in the previous section, however, now, what will 
be interesting is the correlations across the horizon, considering that the different graph 
sectors can transform into each other. Thus the tracing has to be carefully implemented, 
including the correlations across the horizon. Now the density matrix in the previous case 
was derived in a local way, with the tensor product structure giving a diagonal matrix in 
the end. Here the procedure will be quite similar, except that the correlation function / has 
to be implemented in the generalised coherent state. 

IV'gen >= © n l^e„,e«e„, > (37) 

where the above corresponds to the set of edges at the horizon for each inequivalent graph 
considered above. 

Now, the correlation function f{g{o}:9{i},PH) which is introduced to approximate the 
correlated coherent state to implement the 'local apparent horizon equation', now has the 
added complication to ensure that they are independent of the graph F to which the edge 
belongs. As obtained in the previous calculation, if the correlation function at the classical 
level depended on the graph, then it would have to include global properties, such as, for 
example, the total number of edges crossing the horizon, the distribution of the vertices in 
the classical metric, etc, which clearly is not true. 

® n fg{o}^PH,g{i} 

In the classical limit, the derivation of the density matrix will follow the exact procedure 
as above, and the off-diagonal elements for different graphs would be damped in the semi- 
classical limit. Thus in the t — ^ limit the density matrix will comprise only of the diagonal 
elements, but now, each of them proportional to the same function The density ma- 

trix thus obtained is a Z^jjj n-t(2ji + 1) x J2{ji}Y[i{'^ji + 1) matrix. The sum is over all 
possible configurations of the horizon graph. The density matrix has diagonal entries, each 
proportional to |/p hence, from Trp = 1, one obtains the individual diagonal entries to be 
1/ HiC^ji + 1))) • The entropy obtained from this density matrix is then 

S = -Tr(plnp) = ln(^n(2j. + 1)) (39) 
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The combinatorial counting for this is calculated exactly in the next section. This entropy 
is again proportional to area, and correctly gives Bekenstein Hawking entropy. 

A. Combinatorial Counting 

In this section, we consider three enumerative problems. First is the well-understood 
question of determining d, as defined in (33), the number of different ways there are to 
choose a positive integer n and spins ji for 1 <i <n subject to the constraint that 

E(J. + 1/2) = A. (40) 

i 

This constraint is due to the bits of semiclassical areas induced by the edges which have to 
add up to the total area A of any macroscopic surface measured in units of Planck length 
squared. Secondly, we approximate \ Y,jA'^^i{'^3i + 1)); where the sum is taken over all 
assignments of n and ji satisfying (40). We show that this sum is bounded by a constant 
not depending on A. 

Finally, we give an exact calculation of In Y^ji Wii'^ji + 1) which is required for the solution 
of (39). 

For the duration of this section, we write for 2(jj + l/2) = 2jj + l. The requirement that 
the ji be non-negative half-integers means that the are required to be positive integers. 
We also write for 2A. 

The first problem, then, is how many ways to choose n and Oj for 1 < i < n positive 
integers, such that X^i o-i = (These are called compositions of in the combinatorics 
literature.) To solve this, imagine that we have a real line of length A^, marked off into 
equal pieces by A^ — 1 perpendicular lines. If we choose any subset of the lines as positions to 
cut the real line, and then consider (in order) the lengths of the smaller pieces obtained, we 
will obtain a possible assignment of values for where n is the number of pieces. 

Thus, our problem is equivalent to asking how many ways we can choose a subset of the A^ — 1 
lines to cut. Each line is either chosen or not chosen, so the total number of compositions 
oi N is d = 2^-\ 

Now we consider the second problem. 

N 

EMn(«^)) = EMA^)ln(A;), (41) 

at i k=l 
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where /i(/c) is the total number of times k appears as a part, among all the d compositions 
of N. 

Clearly, fi{N) = 1, since there is only one composition in which appears as a part, 
namely N = N. To understand fi{k) for k < N, consider fir{k), where fir{k) is the number 
of times k appears as a summand in one of the d compositions, such that the total of the 
summands appearing before k is exactly r. 

Thus, fio{k) counts the number of compositions of where the first part is k. Clearly, 
by our solution to the first problem, fio{k) = 2^"^^^. Also, ^N-k{k) counts the number of 
compositions of A^ whose final part is fc, and thus, by symmetry, fiN-kik) = 2^"^^^. 

Now consider fir{k) for < r < A^. We see that it counts compositions of A^ which consist 
of a composition of r, followed by a part of size k, followed by a composition of N — k — r. 
Thus, 

firik) = 2'-i2^-'="^"i = 2^-'=-^. (42) 
Summing, we find that for any 1 < k < N, 

N-k 

^(k) = ^ fxrik) = {N-k + 3)2^-'=-2 (43) 

r=0 

Thus, the quantity we are interested in is given by 

(j:k-oiN-k + 3)2^~^-Hnik))+ln{N) / , , \ 

^ ° — — = E '2 (N -k + 3) \nik) + 2"^+i In(A^). 

(44) 

The negative exponential predominates here, with the result that as A^ goes to infinity, 
the sum stays bounded. 

Now we consider the third problem. We want to understand 

n ai + ...+an=N i 

The second sum is over all compositions of A^. To help us understand P/v, let us fix n 
and write 

PN,n= E (45) 

ai+...+a„=N i 

We encode the values PN,n in a generating function: 

oo 

UX) = E PN,nX''. (46) 
Af=0 
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Now we claim that = (x + + 3x^ + Ax'^ + . . .)"". Compare x^ coefficients on both 

sides. The x^ coefficient on the lefthand side is P/v,n by assumption. On the righthand side, 
contributions to will be obtained by choosing a term from each of the n factors, such 
that the exponents of the chosen terms add to A^. Thus, this choice amounts to a choice of a 
composition of A^. Now we see that the coefficients were chosen precisely so that the weight 
corresponding to a given composition is the product of the parts, and thus the coefficient of 
on the righthand side is also PN,n, as desired. 

So /„(x) = (x + 2x^ + 3x'^ + ...)" = . Summing, we find that 



Thus, we have obtained a generating function for Pn- Note that this is essentially the 
same calculation as was carried out in [22]; their partition function Z is our generating 
function /(x) once one sets = x. 

Write 1 — 3x + x^ = (1 — ax)(l — /5x), where 

3 + V5 ^ 3-V5 
2 ' ^ 2 

Expanding in partial fractions, we get 

^ iZ^ - (48) 




Therefore, 



3 + v^\^ /3-v^\^^ 



1 //3 + v^\^ /3-v^\^^ 



x^ (49) 



Thus, the quantity we are interested in, \n{P2A), is given by 

2A 

/ 'J . / R \ \ \ 

(50) 

Thus the above would the the degeneracy of the area of any macroscopic surface with 
total area A. 
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B. Entropy 

Now, after obtaining exact expressions for the various ways of counting the horizon, we 
now speciahse to the determination of entropy. The basic formahsm is the same: we are 
counting the number of ways to induce the horizon area by arbitrary graphs, in the semi- 
classical limit. The horizon area is given by a sum of bits of areas which equal some integer 
given by (jci + where jci is the spin of the edge inducing the area. 

In most of the previous countings of entropy, the spins were all set equal to 1/2, and a 
justification for that appears in [29]. However, here, they would correspond to the choice 
of a very symmetric graph, where all the spins induce the horizon with precisely an area t. 
The entropy counting would thus be trivial and equal 

'5BH = ^ln2 (51) 

Now, this is not exactly the same as the Bekenstein-Hawking entropy as the constant is 
not 1/4. This problem is solved as in [11], by realising that the area apectrum is modified 
by the Immirzi parameter in different quantisation sectors of the theory. Thus the area 
bit is (j + l/2)tj3, and hence the entropy term gets modified now by the requirement that 
the integers have to add up to area divided by the Immirzi parameter. Thus the constant 
= 2A/(3 (in units of Planck length squared). The Immirzi parameter would be 4 In 2, 
though not exactly equal to the value obtained from previous calculations of entropy [11]. 
In fact, what is being identified here is a problem: the lack of universality of the 'Immirzi 
parameter'. Every time one gives a counting of black hole entropy within the formalism 
of loop quantum gravity, one discovers a new value of the same parameter. Though by 
definition, the parameter is an ambiguity in quantisation, these 'theoretical' experiments do 
not actually fix the ambiguity in any way. Firstly, given the level of theoretical knowledge, 
it is difficult to determine the correct entropy calculation. However, the universality of 
Bekenstein-Hawking entropy as being proportional to the area of the horizon is true for all 
countings. This law will hopefully be verifiable experimentally, by actual observations. Note 
also in the above case (51), the counting is exact, and there are no corrections semi-classically 
to the area law. There will be quantum fluctuations which will correct the classical space- 
time, however that need not be proportional to area and would be one order higher in the 
semi-classical parameter t. 
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Now, suppose the number of spins were i = l..n. Then the density matrices would be 
nr=i(2j« + 1) X nr=i(2jj + l) matrices and hence the entropy would be of the form: 



Sbh = EM2j^ + l) 

i 

= ln(n(2j. + l)). 



(52) 
(53) 



In addition there will be the constraint that the spins sum up to the total area I]i(ji + l/2)^ = 
An/lp- For a fixed graph, where the number of edges is fixed, the classical area Ph would 
thus be fixed and hence the spin counting i = l..n also fixed to be a certain number. 
However, what is now different from the all spins equal case, is that area is not equally 
distributed, and hence in some sense one is summing over inequivalent graphs but with the 
same number of edges crossing the horizon. This would correspond to the diagram of the 
previous section. This entropy would correspond to the log of P/v,„ of equation (45). To 
extract the coefficient of PN,n from 

DO / 



U{x) = x"(l - x) 



-2n 



i=0 



\ 



-2n 
i 



(-l)V 



(54) 



Hence, the coefficient of x will be 



N,n 



-2n 
N -n 



-1) 



N-n 



2n{2n + l) {n + N - 1] 

{N-n)\ 



(55) 



Now, the n can be absolutely arbitrary, but that will create a dependence of the entropy on 
this unknown parameter. However, to gauge some understanding of the entropy we take the 
number of spins to be N/2, so as to compare with the case where all the spins were set to 
1/2. 

Clearly, 

P„„_^^(^^±4.fc).4.-^ (56) 



N,N/2 



2 ^' 



N 



2-: 3A^(Ar-l)!f! 



Using the Stirling's formula, r! = r'^e \j2'nr, for large r, and one obtains the following: 

~NI2 



Pn,N,2 - ^{ — ) e Y - 1) 



V27riVv^ 



2^3 



3V2^iVV 



33^/22-^ + 0(1) 

1/2 



(57) 
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The entropy or log of this quantity is then 

Sbh = N{- In 3 — In 2) In H — ln( — ) + further corrections (58) 

Now, = 2^, and thus we are observing a particular value of the Immirzi parameter, 
and a correction to Bekenstein-Hawking [2, 3] entropy which is proportional to log of area, 
already observed. Of course, here, the area spectrum is different, and the counting is exact. 
Note that this entropy is different from the entropy when all the spins were symmetrically 
distributed (51). 



C. Sum over all spins 

While trying to understand the sum over all possible graphs, and the corresponding 
entropy, one takes, as a first case, the tensor sum of density matrices for individual graphs. 
The entropy is then given by (33). The combinatorial counting for d has been done, and the 
correction term is exponentially supressed as given in (44). Thus the entropy is 

Sbh ~ {N - l)\n2 + exp{-{N - l)\n2)\nN (59) 

Note that = 2Ah/Ip/3 and there is indeed a In A term, however, that is multiplied by a 
term which decreases exponentially with area. If one expanded the exponential, then the 
first term would be a log correction term, however, this would be true only for small enough 
areas. 

Now, we take the physically interesting case where density matrix is ab-initio derived 
from the generalised coherent state. This is given by the expression for P2A in equation 
(50). This is surprisingly similar to the one obtained in [10], where the counting was done 
for the number of integer partitions of a given area. The entropy is thus 

lnP2A = 2^1n(^±^) +ln(l - (6.854)-'^- ^K^) -Unb (60) 

Pv 

= 2^ ln(^±^) + (6.854) -2^-/^'^^ + ... (61) 

Clearly, the first term is the Bekenstein-Hawking entropy, but the corrections now are ex- 
ponentially decreasing in area. As area increases, the corrections go to zero. The Immirzi 
parameter obtained here is precisely proportional to the one obtained in [22], which is not 
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surprising as the counting here is for the same object. This new observation regarding 
correction to entropy might have a origin from quantum information theory [30]. The cor- 
rection to the entropy is a completely new result here, though the implications of this are not 
very clear. Classical laws of black hole mechanics are usually supplemented with Hawking 
temperature to get the exact 1/4 factor. But this can have semi-classical corrections and 
considerably modify our interpretations. 

IV. CONCLUSIONS 

We defined generalised coherent states over a tensor sum of Hilbert spaces labelled by 
minimal graphs and their corresponding classical labels. We derived the density matrix from 
the corresponding generalised coherent state and showed that the entropy is proportional 
to area. The corrections to entropy are exponentially decreasing with area. These coherent 
states are still defined in the kinematic Hilbert space. With the new work on the Hamiltonian 
constraint, we hope to be able to contribute to the derivation of the coherent state in the 
constraint physical Hilbert space. However, since in the limit t ^ 0, the results predicted 
by those physical coherent states have to agree with the classical results, we do not expect 
much deviation from the conclusions of our paper. 
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